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Abstract 

Cascade is a full hadron level Monte Carlo event generator for ep, 7p and pp processes, 
which uses the CCFM evolution equation for the initial state cascade in a backward 
evolution approach supplemented with off - shell matrix elements for the hard scattering. 
A detailed program description is given, with emphasis on parameters the user wants to 
change and common block variables which completely specify the generated events. 

1 Tabular Summary 



program name 
version 

date of latest version 
author 
program size 
input files needed 
computer types 
operating systems 
applicability 



hard sub-processes included 



QCD cascade 

fragmentation model 
other programs called 

availability 



Cascade 
1.00/01 
July 2001 

Hannes Jung (Hannes.Jung@desy.de) 
~ 8000 lines of code 
ccfm.dat (kms.dat, kmr.dat) 

any with standard Fortran 77, tested on SGI, HP-UX, SUN, PC 
Unix, Linux 

Deep Inelastic ep Scattering 
Photo-production in ep Scattering 
7P Scattering 

heavy quark production in pp and pp Scattering 

1*9* qq, QQ 
9*9* qq. QQ 

initial state parton shower according to CCFM 
final state parton shower with angular ordering 
LUND string 
Pythia 6.1 
Bases/Spring 5.1 

http : //www. quark. lu. se/~hannes/ cascade/ 
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Figure 1: Kinematic variables for multi-gluon emission. The t-channel gluon four-vectors are 
given by ki and the gluons emitted in the initial state cascade have four-vectors pi. The maxi- 
mum angle for any emission is obtained from the quark box, as indicated with H. 



The formulation of the CCFM [ ^ ^ ^, ^ parton evolution for the implementation into a full 
hadron level Monte Carlo program is described in detail in [ 0. Here only the main results 
are summarized and discussed for the case of leptoproduction. 

Figure |l| shows the pattern of QCD initial-state radiation in a small-x DIS event, together 
with labels for the kinematics. According to the CCFM evolution equation, the emission of 
partons during the initial cascade is only allowed in an angular-ordered region of phase space. 
The maximum allowed angle S for any gluon emission sets the scale q for the evolution and is 
defined by the hard scattering quark box, which connects the exchanged gluon to the virtual 
photon. In terms of Sudakov variables the quark pair momentum is written as: 



where Pp and Pe are the proton and electron momenta, respectively and Qt is the transverse 
momentum of the quark pair and Y is its light-cone momentum fraction. The scale q is related 



for a backward evolution approach adopted in the Monte Carlo generator Cascade [ ^: 



2 The CCFM evolution equation 



Pg+Pq = Y{pp + Epe) +Q, 



(1) 




_2 d xA{x, kt, q) 



dz 



dxp P{z,q/z,kt) 
2tx As{q,Qo) 



x'A{x', k[, q/z) 



(2) 



dq^ A,(g,Qo) 
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where A{x, kt, q) is the unintegrated gluon density, depending on x, kt and the evolution variable 

q. The splitting variable is z = x/x' and kt = {\ — z)/zq + kt, where the vector q is at an 
azimuthal angle 0. The Sudakov form factor is given by: 

A,(,, Q„) = exp I dz j (3) 

with ots = = For inclusive quantities at leading-logarithmic order the Sudakov form 
factor cancels against the 1/(1 — z) coUinear singularity of the splitting function. 

The splitting function Pg{zi, qi, kti) for branching i is given by: 

P,{z,. g., ktd = "^^f^"'^^'^ + ^A„.(.„ ql kl) (4) 
where the non-Sudakov form factor A^^ is defined as: 

log = /' ^ / %<^{kt. - g)e(g - z'g,.) (5) 

Jo z J q 

The splitting function P(z, g, kt) in eq.(^ contains only the singular parts in z and {\ — z). 
The finite terms in the splitting function are neglected, since they are not obtained in CCFM 
at the leading infrared accuracy [ ^, p. 72]. As soon as they become available they can be easily 
implemented into the Monte Carlo generator. 



3 Backward evolution: CCFM and CASCADE 

The idea of a backward evolution [ 0, |^ is to first generate the hard scattering process with 
the initial parton momenta distributed according to the parton distribution functions. This 
involves in general only a fixed number of degrees of freedom, and the hard scattering process 
can be generated quite efficiently. The initial state cascade is generated by going backwards 
from the hard scattering process towards the beam particles. 

According to the CCFM equation the probability of finding a gluon in the proton depends 
on three variables, the momentum fraction x, the transverse momentum squared k\ of the 
exchanged gluons and the scale q = x„a/sS, which is related to the maximum angle allowed 
for any emission H. To solve eq.( ||) the unintegrated gluon distribution A{x,k,q) has to be 
determined beforehand. 

Given this distribution, the generation of a full hadronic event is separated into three steps, 
as implemented in the hadron-level Monte Carlo program Cascade: 

• The hard scattering process is generated, 

a = / dktdXgA{xg, h, q)cr{Y9* X) , (6) 
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with X being gg, QQ or Jjiljg states. The hard cross section is calculated using the off- 
shell matrix elements given in [P, p. 178 ff] for qq and QQ or for the case of 7(7* — > J/ipg 
rO| . The gluon momentum is given in Sudakov representation: 



m 



k = XgPp + XgPe + kt^ XgPp + /Cf . (7) 

where the last expression comes from the high energy approximation {xg <^ 1), which 
then gives —k"^ ~ /cj . 

The initial state cascade is generated according to CCFM in a backward evolution ap- 
proach. 

The hadronization is performed using the Lund string fragmentation implemented in 



Pythia /Jetset[ ITT 



The backward evolution there faces one difficulty: The gluon virtuality enters in the hard 
scattering process and also influences the kinematics of the produced quarks and therefore the 
maximum angle allowed for any further emission in the initial state cascade. This virtuality is 
only known after the whole cascade has been generated, since it depends on the history of the 
gluon evolution (as Xg in eq.( |^) may not be neglected for exact kinematics). In the evolution 
equations itself it does not enter, since there only the longitudinal energy fractions z and the 
transverse momenta are involved. This problem can only approximately be overcome by using 
fc^ = k^/{l — Xg) for the virtuality which is correct in the case of no further gluon emission in 
the initial state. 

The Monte Carlo program Cascade can be used to generate unweighted full hadron-level 
events, including initial-state parton evolution according to the CCFM equation and the off- 
shell matrix elements for the hard scattering process. It is applicable for pp, photo-production as 
well as for deep inelastic scattering. The typical time needed to generate one event is ~ 0.03 sec 
on a Pentium II (266 MHz), which is similar to the time needed by standard Monte Carlo event 
generators such as Lepto [ O or Pythia [ |n|. 



3.1 The unintegrated gluon density 

The CCFM unintegrated gluon density xA{x, kt, q) is obtained from a forward evolution pro- 
cedure as implemented in Smallx [ |T3|, [l^ by a fit to the measured structure function F2 as 
described in [ |, 0. From the initial gluon distribution a set of values x and kt are obtained by 
evolving up to a given scale log q using the forward evolution procedure. This is repeated 10* 
times thus obtaining a distribution of the unintegrated gluon density xA{x, kt , q) for the slice 
of phase space with a given q {q > qtn)- To obtain a distribution in logg, the above procedure is 
repeated from the beginning 50 times for the different grid points in log g up to g = 1800 GeV. 
Due to the complicated structure of the CCFM equation, no attempt is made to parameterize 
the unintegrated gluon density. Instead, the gluon density is calculated on a grid in log x, log kt 
and log g with 50 x 50 x 50 points and a linear interpolation is used to obtain the gluon density 
for values in between the grid points. The data file (ccfm.dat) containing the 50 x 50 x 50 
grid points is read in at the beginning of the program. 

The parameter IGLU offers the possibility to compare the CCFM unintegrated gluon density 
xA{x,kt,q) (IGLU=1) with other unintegrated gluon densities published: a simple numerical 
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derivative of a standard integrated gluon density '^^^^^ ^ taken from [ ^] (IGLU=2), the one in 
the approach of Bliimlein [ |16| and coded by [ 0, |18[ (IGLU=3), the unintegrated gluon density 
of KMSQ [ |ll (IGLU=4, stored in kms .dat), the one of the saturation model by [ |^ (IGLU=5) 
and the one of KMRQ [ ^ (IGLU=6, stored in kmr.dat). 

Except for the default CCFM unintegrated gluon density with (IGLU=1), no initial state 
parton shower can be generated, because the angular variable, essential for angular ordering 
in the initial state cascade, is not present. However, the transverse momenta of the incoming 
partons are properly treated. 



3.2 as and the choice of scales 



The strong coupling is calculated via the Pythia [ subroutine PYALPS. Maximal and 
minimal number of flavors used in are set by MSTU(113) ,MSTU(114) , Kqcd =PARU(112) 
with respect to MSTU(112) flavors and stored in the Pythia common COMMON/PYDATl/. In the 
initial state cascade according to CCFM, the transverse momenta of the t-channel gluons are 
allowed to perform a random walk for small z values and kt can become very small. In the 
1/ z part of the splitting function we use fi = kt a.s the scale in as{fj.) and in the 1/(1 — z) part 
/i = (1 — z)qt is used. In addition we require fi > Qo = 1.4 GeV, resulting in ^^(/i > Qq) < 0.6. 
The parameter Qq = 1.4 GeV was determined from the requirement to describe the structure 
function F2 as described in [ ^ Q . 

The scale /i which is used in as in the hard scattering matrix element can be changed with 
the parameter IQ2, the default choice is //^ = pf. 



3.3 Quark masses 

The quark mass for light quarks {u,d,s) is fixed to ruq = 0.140 GeV. This, together with the 
treatment of at small scales /i, gives also a reasonable total cross-section for photo-production 
at HERA energies. The masses for heavy quarks are given by the Jetset / Pythia defaults 
{nic = 1.5 GeV, mb = 4.8 GeV) and can be changed according to the Pythia prescription. 



3.4 Leptoproduction 

Cascade can be used to simulate leptoproduction events over the whole range. By fixing 
the light quark masses to = 0.140 GeV and for small /i, the hard scattering matrix 
element remains finite over the full phase space. The total cross section is simulated by se- 
lecting IPR0=10 and NFLAV=4(5). With IPR0=10 light quarks {u,d,s) are selected and with 
NFLAV>3 the program automatically includes heavy flavor production via the process IPR0=11 
and IHFLA=4 up to IHFLA=NFLAV. 

Heavy flavor production can be generated separately via IPR0=11. The value of IHFLA 
determines the heavy flavor to be generated. 



The matrix element for 'jg* — > J/ipS calculated in [[T^ is available for quasi-real 7's via the 
process IPR0=2. 



^A. Stasto kindly provided the program. 
^M. Kimber kindly provided the program. 
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3.5 Photoproduction 



Cascade can be used to simulate real photoproduction events by using KE=22. The same 
options as for leptoproduction are available. 



3.6 Hadroproduction 

Cascade can be used to simulate heavy quark production in pp or pp collisions (IPR0=14 for 
heavy flavor production, and IHFLA=4(5) for charm (bottom) quarks). The flavor code (KE) for 
beam 1 can be chosen as KE=2212 for proton or KE=-2212 for anti - proton, beam 2 is always 
a proton. 



3.7 Initial and final state parton shower 

Initial state parton showers are generated in a backward evolution approach described in detail 
in [ ^, . The initial state parton shower consists only of gluon branchings and is generated in 
an angular ordered region in the laboratory frame. In the present version, the gluons emitted 
during the branchings are treated on-mass shell, and no further time - like branching occurs. 

All parameters (like the scale fi in as, the coUinear cut-off Qq) for the initial state cascade 
are fixed from the determination of the unintegrated gluon density. The transverse momenta 
of the partons which enter the hard scattering matrix element are already generated in the 
beginning and are not changed during the whole initial and final state parton showering (in 
contrast to standard DGLAP type parton shower Monte Carlo generators like Pythia [ |1^ 
or Lepto [IH). 

The final state parton shower uses the parton shower routine PYSHOW of Pythia with the 
scale QMAX=(mi ± + m2 ±)/2, with mi(2) _l being the transverse mass of the hard parton 1(2). 



3.8 Remnant treatment 

The proton remnant is built in subroutine CAREMN, which is a slightly modified version of the 
Pythia / Lepto subroutine PYREMN. In the present version no intrinsic transverse momentum, 
in addition to the transverse momentum from the initial state cascade, is included. 



4 Description of the program components 

In Cascade all variables are declared Double Precision. The Lund string model is used for 
hadronization as implemented in Pythia 6.1 [ 11]. The final state QCD radiation is performed 
via PYSHOW from Pythia 6.1. The treatment of the proton remnant follows very closely the 
ones in Lepto[ |1^ for the leptoproduction case and the one in Pythia 6.1 for the proton - 



proton case. However slight modifications were needed to adapt to the cascade treatment here. 

The unintegrated gluon density is stored on data files (ccfm.dat,kms.dat,kmr.dat), and 
is read in at the beginning of the program. 

The program has to be loaded together with Pythia 6, to ensure that the double precision 
code of Jetset is loaded. 
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4.1 Random number generator 



Since the variables are declared in double precision, also a double precision random number 
generator has to be used to avoid any bias. The function DCASRN gives a single random number, 
the function DCASRNV returns an array of length LEN of random numbers. The default random 
number generator is RM48 (called in DCASRN and DCASRNV) from Cernlib. The user can change 
this to any preferred Double Precision random number generator. 



4.2 Integration and event generation 



The integration of the total cross section and the generation of unweighted events is performed 
with the help of Bases/Spring [ ^|. 



4.3 Subroutines and functions 



The source code of Cascade and this manual can be found under: 
http : //www . quark . lu . se/~hannes/ cascade/ 

CAMAIN main program. 

CASINI to initializes the program. 

CASCADE to perform integration of the cross section. This routine has to be called 

before event generation can start. 
CAEND to print the cross section and the number of events. 

CAUNIGLU(KF , X , KT , P , XPQ) to extract the unintegrated gluon density xA{x, kf, q) for a proton 
with KF=2212, as a function of x =X, =KT and q =P. The gluon density is 
returned in XPQ(O), where XPQ is an array with XPQ (-6: 6). 
to perform the event generation, 
to give with =RQ. 

to give initial particle and parton momenta, 
to call routines for selected processes: XSECl. 
to cut on pt for 2-^2 process in integration and event generation, 
matrix element for 7*5'* qq and g*g* — qq including masses, q can be 
light or heavy quarks, 
four- vector dot product of A and B. 

four-vector dot product of vectors I and J in PYJETS common, 
to generate momenta of final partons in a 2 ^ 2 subprocess according to 
phase space 

P_SEMIH to generate kinematics and the event record for ep, 7p and pp processes. 

CAREMN(IPU1,IPU2) to generate the beam remnants. Copied from LEPTO 6.1 [ [12| and 

updated for the use in CASCADE. 
CASPLl (KF , KPA , KFSP , KFCH) to give the spectator KFSP and KFCH partons when a parton KPA 



EVENT 
ALPHAS (RQ) 
PARTI 
FXNl 

CUTG(IPRO) 
MEOFFSH 

DOT(A,B) 

DOTKI.J) 

PHASE 



is removed from particle KF. Copied from LEPTO 6.1 [ |T2| and updated for 
the use Cascade. 

CAPS to generate color flow for all processes and prepare for initial and final state 

parton showers. 
CASCPS (IPUl , IPU2) to generate initial state radiation. 

GADAP Gaussian integration routine for 1-dim and 2-dim integration. Copied from 



LEPTO 6.1 
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4.4 Parameter switches 

BASES/SPRING Integration procedure. 
NCAL : (D:=20000) Nr of calls per iteration for bases. 

ACCl : (D:=l) relative precision (in %) for grid optimization. 

ACC2: (D:=0.5) relative precision (in %) for integration. 



4.4.1 Parameters for kinematics 



PLEPIN: (D:=— 30) momentum p [GeV/c] of incoming electron (/INPU/). 

PIN: (D:=820 momentum p [GeV/c] of incoming proton (/INPU/). 

QMI: (D: = 5.0) (/VALUES/) minimum Q'^ to be generated. 

QMA: (D: = 10^) (/VALUES/) maximum to be generated. 

YMI : (D: = 0.0) (/VALUES/) minimum y to be generated. 

YMA: (D: = 1.0) (/VALUES/) maximum y to be generated. 

THEMA,THEMI (D: THEMA = ISO., THEMI = 0) maximum and minimum scattering angle 6 of 

the electron (/CAELEC/). 

NFLAV (D: = 5) number of active flavors, can be set by user (/CALUCO/). 



4.4.2 Parameters for hard subprocess selection 

IPRO: (D: = 10) (/CAPARl/) selects hard subprocess to be generated. 

-2: -fg* ^ J/ijg 

—10: 7**7* qq for light quarks. 

=11: Yg* ~^ QQ for heavy quarks. 
PT2CUT(IPR0) : (D=0.0) minimum for process IPRO (/CAPTCUT/). 



4.4.3 Parameters for parton shower and fragmentation 

NFRAG: (D: = 1) switch for fragmentation(/CAINPU/). 

= off 

= 1 on 

IFPS: (D: = 3) switch parton shower (/CAINPU/). 

= off 

= 1 initial state 
= 2 final state 
= 3 initial and final state 
ICCFM: (D: =1) 

=1 CCFM evolution. 



4.4.4 Parameters for structure functions and scales 

IRUNAEM: (D: = 0) (/CAPARl/) select running of aemiQ'^)- 

=0: no running of aem{Q^) 

=1: running of Q!em(Q^) 
IRUNA: (D: = 1) switch for running a^. 

=0: fixed cts = 0.3 
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IQ2: 



IGLU: 



=1: running a<((/i^) 

(D: = 3) select scale /i^ for Q;s(/i^). 

=1: = 4 ■ (use only for heavy quarks!) 

=2: fi"^ = s (use only for heavy quarks!) 

=3: jj,"^ = 4 ■ m'^ + p]_ 

=4: fi' = Q' 

=5: ^^ = Q^+pI 

(D: = 1) select unintegrated gluon density/GLUON/). 

Note for IGLU> 1, all initial state parton showers are switched off. 



CCFM unintegrated gluon (ccfm.dat) 

dxg{x,Q'^) 
dQ2 • 



derivative of GRV 
approach of Bliimlein 
KMS [|191 (kms.dat). 
saturation model [ ^ 
KMR \ M (kmr . dat) 



m 



4.4.5 Accessing information 



AVGI integrated cross section (/CAEFFIC/). 

SD standard deviation of integrated cross section (/CAEFFIC/). 

SSS squared center of mass energy s (/CAPARTON/). 

PBEAM energy momentum vector of beam particles (/CABEAM/). 

KBEAM flavor code of beam particles (/CABEAM/). 

Q2 in leptoproduction: actual of 7 (/CAPAR4/). 

YY negative light-cone momentum fraction of parton 1 (7*, g*) (/CASGKI/). 

YY_BAR positive light-cone momentum fraction parton 1 (7*, g*) (/CASGKI/). 

XG positive light-cone momentum fraction of parton 2 (g*) (/CASGKI/). 

XG_BAR negative light-cone momentum fraction of parton 2 (g*) (/CASGKI/). 

KT2_1,KT2_2 transverse momenta squared k^^^^ [GeV^] of partons 1(2) which enter to the 
matrix element. 

YMAX,YMIN actual upper and lower limits for y =YY (/CAPAR5/). 

Q2MAX , Q2MIN actual upper and lower limits for Q'^ (corresponding to KT2_l) of 7 (/CAPAR5/). 

XMAX,XMIN upper and lower limits for x (/CAPAR5/). 

AM (18) vector of masses of final state particles of hard interaction (/CAPAR3/). 

SHAT invariant mass s [GeV^] of hard subprocess (/CAPAR5/). 

NIA1,NIA2 position of partons in hard interaction in PYJETS event record (/CAHARD/). 

NF1,NF2 first and last position final partons/particles of hard interaction in PYJETS 
(/CAHARD/). 

Q2Q hard scattering scale /i^ used in and structure functions (/CAPAR4/). 

ALPHS actual (/CAPAR2/). 

ALPH aem (/CAPAR2/). 

NIN number of trials for event generation (/CAEFFIC/). 

NOUT number of successful generated events (/CAEFFIC/). 
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4.5 List of COMMON blocks 



C0MM0N/CABEAM/PBEAM(2,5) ,KBEAM(2,5) ,KINT(2,5) 
COMMON/CAHARD/ NIAl , NIA2 , NIR2 , NFl , NF2 
COMMON/CAHFLAV/ IHFLA 

CQMMON/CAINPU/PLEPIN , PPIN , NFRAG , ILEPTO , IFPS , IHF , INTER , ISEMIH 
COMMON/CALUCO/KE , KP , KEB , KPH , KGL , KPA , NFLAV 
COMMON/CAEFFIC/AVGI , SD , NIN , NOUT 
COMMON/CAELEC/THEMA , THEMI 
COMMON/CAGLUON/IGLU 

COMMON/CAPARl/IPRO , IRUNA , IQ2 , IRUNAEM 

COMMON/CAP AR2/ALPHS , PI , ALPH , IWEI 
C0MM0N/CAPAR3/AM(18) , PCM (4, 18) 
COMMON/CAP AR4/Q2 , Q2Q 

C0MM0N/CAPAR5/SHAT , YMAX , YMIN , Q2MAX , Q2MIN , XMAX , XMIN 
CQMM0N/CAPAR6/LST(30) ,IRES(2) 
COMMON/CAP ARTON/SSS , CM (4) , DBCMS (4) 
C0MM0N/CAPTCUT/PT2CUT (20) 

COMMON/CASKIN/YY , YY.BAR , XG , XG.BAR , KT2_ 1 , KT2_2 , PT2H , SHH 
COMMON/VALUES/QMI , YMI , QMA , YMA 



5 Example Program 



PROGRAM CASMAIN 
Implicit None 
Integer N1,N2 
REAL PLEPIN,PP1N 

INTEGER KE , KP , KEB , KPH , KGL , KPA , NFRAG , ILEPTO , IFPS , IHF 
INTEGER INTER, ISEMIH 

INTEGER NIAl , NIRl , NIA2 , NIR2 , NFl , NF2 , NET , NFLAV 
CGMMON/CALUCG/KE , KP , KEB , KPH , KGL , KPA , NFLAV 

COMMON/CAINPU/PLEPIN , PPIN , NFRAG , ILEPTO , IFPS , IHF , INTER, ISEMIH 
COMMON/CAHARD/ NIAl , NIA2 , NIR2 , NFl , NF2 
INTEGER IHFLA 
COMMON/CAHFLAV/ IHFLA 

DOUBLE PRECISION THEMA , THEMI , PT2CUT 
INTEGER IRUNA , IQ2 , IRUNAEM 
INTEGER IPRO 

COMMON/CAPARl/IPRO , IRUNA , IQ2 , IRUNAEM 
COMMON/CAELEC/ THEMA, THEMI 
C0MM0N/CAPTCUT/PT2CUT (20) 
REAL ULALPS,ULALEM 
EXTERNAL ULALPS , ULALEM 
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DOUBLE PRECISION QMI , YMI , QMA , YMA 
COMMON/VALUES/QMI , YMI , QMA , YMA 

Integer Iglu 
Common/CAGLUON/ Iglu 

Integer I SEED, I 



I SEED = 124567 

nl=0 

n2=0 

C initialize random number generator 
CALL RM48IN(ISEED,N1,N2) 

C initialize PYTHIA 6 parameters 
CALL GPYINI 

C initialize CASCADE parameters 
CALL CASINI 



C Select parton shower (IPS=1 initial, =2 final, 3 initial+f inal PS ) 

IFPS = 3 
C scale for alpha_s 

C IQ2 =1 mu"2 = m_q~2 (m_q = light quark or heavy quark depending on IPRO) 
C IQ2 =2 mu''2 = shat 

C IQ2 =3 mu''2 = m_q"2 + pt''2 (m_q = light quark or heavy quark depending on IPRO) 

C IQ2 =4 mu"2 = q~2 (q~2 of virtual photon) 
C IQ2 =5 mu"2 = q~2 + pt"2 (q"2 of virtual photon) 
IQ2=3 

C select process (IPR0=10 for light quarks, IPR0=11 for heavy quarks) 
IPRO= 10 

C total number of flavors involved 
NFLAV = 4 

C select unintegrated gluon density (D=l) 
Iglu = 1 

C minimum Q''2 of electron to be generated 

QMI = 0.5d0 
C maximum Q"2 of electron to be generated 

QMA = 10D8 
C minimum y of electron to be generated 

YMI=0 . OdO 

C minimum y of electron to be generated 
YMA=1.0dO 

C maximum theta angle of scattered electron 

THEMA = 180. ODO 
C minimum theta angle of scattered electron 

THEMI = O.ODO 
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C momentum of beam 1 (electron, proton, antiproton) 
PLEPIN =-27.5 

C Lund flavor code for beam 1 (electron=ll,photon=22,proton=2212,antiproton=-2212) 
KE=11 

C momentum of beam 2 (proton) 

PPIN = 820. 
C perform fragmentation NFRAG=0/1 

NFRAG = 1 

c for IPRO = 11 which flavor is produced 
IHFLA = 4 

c 

c Start integration of x-section 
c 

CALL CASCADE 

c 

c Print out result of integration of x-section 
c 

CALL CAEND(l) 

c 

c Start event loop 
c 

Do 1=1,100 
c generate an event 
CALL EVENT 
Enddo 

c 

c Print out of generated events summary 
c 

CALL CAEND(20) 

STOP 
END 
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